Microscopic models of quantum jump super-operators 
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We discuss the quantum jump operation in an open system, and show that jump super-operators 
related to a system under measurement can be derived from the interaction of that system with a 
quantum measurement apparatus. We give two examples for the interaction of a monochromatic 
electromagnetic field in a cavity (the system) with 2-level atoms and with a harmonic oscillator (rep- 
resenting two different kinds of detectors). We show that derived quantum jump super-operators 
have 'nonlinear' form J p = 7 diag [_F(n)apa^i ? (ri)] , where the concrete form of the function F{h) 
depends on assumptions made about the interaction between the system and the detector. Un- 
der certain conditions the asymptotical power-law dependence F(h) = (n + l)~ f3 is obtained. A 
continuous transition to the standard Srinivas-Davies form of the quantum jump super-operator 
(corresponding to j3 = 0) is shown. 

PACS numbers: 42.50.Lc, 03.65.Ta, 03.65.Yz 



I. INTRODUCTION 



In the theory of continuous photodetection and con- 
tinuous measurements the (one-count) quantum jump 
super-operatorfQJS) is an essential part of the formalism 
0,0, S, S, H @, 0,E 0, G3, GH, El , since it accounts for the 
loss of one photon from the electromagnetic field (EM) 
and corresponding photoelectron detection and count- 
ing within the measurement apparatus (MA). One of the 
main equations in this theory is the evolution equation of 
the field's density operator p t , or master equation, which 
reads in the simplest variant as 
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where Hq is the EM field Hamiltonian, 7 is the field- 
MA coupling constant and O is some lowering operator, 
representing the loss of a single photon from the field 
to the environment, that may be detected and counted 
by a duly constructed experimental setup. Defining the 
effective non-hermitian Hamiltonian as [ID, [U, [l5|, [l6| 



H, 
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Eq. ([T]) can be written as (we set here % = 1) 
^jr = -i (HeffPt - PtHl f f ) + jOp t Oi 



(2) 



(3) 



whose formal solution is (see, for example, [H H7j) 
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where 
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ffPo 



PoHl 
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po being the density operator for the field state at t = 0. 
The no-count super-operator exp [L (tk — ifc-i)] evolves 
the initial state po from time tk-i to the latter time tk 
without taking out any photon from the field, it rep- 
resents the field monitoring by a MA. The QJS J» = 
7O • O* is an operation which takes out instantaneously 
one photon from the field. Actually, Tr[Jpo] is the rate 
of photodetection [l8[ . 

The explicit form of the QJS is not predetermined. In 
the phenomenological photon counting theory developed 
by Srinivas and Davies [3 the QJS was introduced ad 
hoc as 



Jsd* = JSDa* a f . 



(5) 



Later, Ben-Aryeh and Brif [l9| and Oliveira et al. [2C 
considered QJS of the form 



Je* = IeE- • E + , 



where 



E_ = (a f a+ 1) 



-1/2, 



and E+ = EL 



(0) 



(7) 



are the exponential phase operators of Susskind and Gl- 
ogower [2lL |22| . These "non-linear" operators allow to 
remove some inconsistencies of the SD theory noticed by 
its authors. 

However, the QJS © was introduced in 0, [2(| also 
ad hoc. Therefore it is desirable to have not only a phe- 
nomenological theory, but also some microscopic models, 
which could justify the phenomenological schemes. The 
simplest example of such a model was considered for the 
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first time in (23[, where the QJS of Srinivas and Davies 
was derived under the assumption of highly efficient de- 
tection. The two fundamental assumptions of that model 
were: (a) infinitesimally small interaction time between 
the field and the MA, and (b) the presence of only few 
photons in the field mode. Only under these conditions 
one can use a simple perturbative approach and arrive 
at the mathematical expression for the QJS, which is in- 
dependent of the details of interaction between the MA 
and the EM. 

If the conditions (a) or (b) are not fulfilled, the QJS 
should depend on many factors, such as, for example, 
the kind of interaction between the field and MA, the 
interaction strength and the time T of the interaction. 
Moreover, it should be emphasized, that the instant tj 
at which the quantum jump occurs cannot be determined 
exactly — it can happen randomly at any moment within 
T. Making different assumptions concerning the moment 
of 'quantum jump', one can obtain different formal ex- 
pressions for the QJS. In [24j we have proposed a simple 
heuristic model for obtaining the 'non-linear' QJS of the 
form 

J» = 7f(a t o)o»a t F(a , a). (8) 

In this connection, the aim of the present paper is to pro- 
vide a more rigorous derivation of QJS's, using a more 
sophisticated model that takes into account dissipation 
effects due to the 'macroscopic part' of the MA. Our 
approach is based on the hypothesis that the transition 
probability must be averaged over the interaction time T, 
during which a photon can be gobbled by the detector at 
any time in the interval (0, T). Considering two different 
models of MA's: a 2- level atom and a harmonic oscil- 
lator interacting with a single-mode EM field, we shall 
demonstrate that different kinds of interaction result in 
quite different QJS's. 

The plan of the paper is as follows. In Sec. II we de- 
rive the QJS using the modified Jaynes-Cummings model 
(with account of damping due to the spontaneous decay 
of the excited state) and calculating the time average of 
the transition operator. In Sec. Ill we apply the same 
scheme to the model of two coupled oscillators, showing 
explicitly how the variation of the relative strength of 
coupling constants results in the change of the function 
F{a)a) in Eq. flSJ. Sec. IV contains a summary and 
conclusions. 



II. MODEL OF TWO-LEVEL ATOM 
DETECTOR 

Let us consider first the model, which is a straightfor- 
ward generalization of the one studied in [231 ] . The role of 
the 'system' is played by a single mode of the electromag- 
netic field, while the 'detector' of the MA (sub-system 
constituting the MA that actually interacts with the EM 
field) consists of a single two-level atom. The Hamilto- 
nian for the total system is chosen in the standard form 



of the Jaynes-Cummings model [25[ 

Ho = -ujQao + Luh + gacr+ + g*a)&- , (9) 

where the Pauli pseudo-spin operators Co and a± cor- 
respond to the atom (<r+ = |e)(p|, cr_ = \g)(e\ and 
(Jo = |e)(e| — |g)(.g|) and one considers that there were 
chosen two levels of the atom (the ground state \g) with 
frequency u) g and the excited state |e) with frequency 
co e = oj g + uio); a, and h = a^a are the lowering, ris- 
ing and number operators, respectively, of the EM field. 
Since the coupling between the field and the atom is 
weak, we assume that to 3> \g\. Until now, the detec- 
tor can absorb and emit photons back into the EM field, 
since the detector is not coupled to some macroscopic 
device that irreversibly absorbs the photons. 

Therefore, we have to take into consideration that the 
detector is coupled to the 'macroscopic part' (MP) of the 
MA (e.g., phototube and associated electronics). Hence 
the detector suffers dissipative effects responsible for the 
spontaneous decay of the excited level of the detector (in 
this case of the atom). And it is precisely this physical 
process that represents a photodetection - the excited 
level of the detector decays, emitting a photoelectron into 
the MP of the MA, which is amplified by appropriate 
electronics and is seen as a macroscopic electrical current 
inside the MP of MA. We can take into account this 
dissipation effects by describing the whole photodetection 
process, including the spontaneous decay, by the master 
equation 

^ + i (HefffH ~ PtHl ff ) = 2Xa^ Pt a+, (10) 

which is the special case of Eq. ((3|), where O = 
O* = er_|_, H e ff = Hq — zA(T + (7_, and 2A is the coupling 
of the excited level of the atom (detector) to the MP of 
the MA (here we make a reasonable assumption that A 
has the same order of magnitude as \g\). The 'sink' term 

R» = 2A(7_ • cr+ (11) 

represents the \e) — > \g) transition within the detector 
(the atomic decay process in this case). If A = 0, then the 
detector interacts with the EM field, but photoelectrons 
are not emitted (thus no counts happen), because the 
absorbed photons are emitted back to the field and then 
reabsorbed at a later time, periodically, analogously to 
the Rabi oscillations. 

In the following, we shall use the quantum trajectories 
approach [l|. The effective Hamiltonian © becomes 

H eff = H - iXa + a_ = - (uj - iX) a 

+um + gao + + g*a? <J- — iX/2 (12) 

(where we have used <j+<j- = (1 + <xo)/2) and the evo- 
lution of the system between two spontaneous decays is 
given by the no-count super-operator 

V tPo = U(t)p U\t), U(t) = exp(-iH eff t). (13) 
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After a standard algebraic manipulation [25j, [26J we ob- 
tain the following explicit form of the non-unitary evolu- 
tion operator U(t): 



U(t) 



exp [— iuj (o~q/2 + n) t] 

Cn+l(t) - i-r-rSn+lit) 

\g\ 



(l + o-o) 



-i T1 Sf i+1 (t)aa + - i-j-rd* Sa+i(t)<T- 

\g\ \g\ 



Cnlt) + l — Sn(t) 

\9\ 



(14) 



where 

C fl (t) = cos (\g\B A t) 7 S A (t) =sm(\g\B A t)/B A , (15) 



B A 



^Jh+(S/\g\f, 5 = \ {w - w - iX) (16) 



(note that parameter 5 is complex and n is an operator). 

Assuming that the field state is po = p F ®\g){g\ at time 
t = or, analogously, the last photoemission occurred 
at t = 0, the probability that the next photoelectron 
emission will occur within the time interval [t, t + At) is 
given by UMM 



P{t) = Tr F _z3 [RDtPo] At, 



(17) 



(the subscripts F and D are a reminder that the trace 
operation is on field and detector spaces, respectively) 
where At is the time resolution of the MA. Tracing out 
first over the detector variables, the probability density 
for the next photoemission to occur at time t will be [27j 



P(*) 



lim 



P(t) 



Ti F [S(t)p F ] , 



(18) 



At^O At 

where the time-dependent transition super- operator 

s(f)« = 2Ar(t).r t (i), (19) 

acting on the EM field, stands for the photoelectron emis- 
sion into the MP of the MA (i.e., the actual photodetec- 
tion) . Once again, the probability for detecting a photo- 
electron in [t,t+ At) is P{t) = Tr F [E(t)p F ] At (now on 
we omit the subscript and write p F = p for the field oper- 
ator). In Eq. (|19[) T(t) is the time-dependent transition 
operator 



T(t) = (e\U(t)\g), 



(20) 



that takes out a single photon from the field state. Sub- 
stituting Eq. (Q3| into Eq. (|2"U| we can write T(t) as 

T(f) = -iA ex p f-Xt/2 - iunt) S A+1 (t)a, (21) 

\g\ 



so the time-dependent transition super-operator (| 19(1 be- 
comes 

E(t)p = 2Xe- xt e-^ fd S fl+1 (t)apa^St +1 (t)e^ M . (22) 



In the resonant case, wo = w, we have 



B,, 



X 



X = X/(2\g\)- 



(23) 



If the interaction time At is small, and the number of 
photons in the field is not very high, in the sense that 
the condition 



\g\AtVr 



1 < 1 



(24) 



is fulfilled for all eigenvalues of n, for which the prob- 
abilities p„ = (n\p\n) are important, then one can re- 
place the operator sin (Bn+i\g\At) in Eq. (|15j) simply by 
Bn+i\g\At and arrive at the QJS 



J p = e 



— iuhAt 



2\{\g\Atf apo) 



(25) 



which has almost the Srinivas-Davies form ([S]), with the 
coupling constant 



1S D = 2\(\g\Atf 



(26) 



Taking 2A = (At) -1 we obtain the same coupling con- 
stant 75 u as in [23| . but this assumption is not the 
only possible. Note that super-operator (|2"B")) contains 
the factors exp(±io;nAt), which can be essentially dif- 
ferent from the unit operator even under condition (|24[) . 
for two reasons: (1) the condition \g\At <C 1 does not 
imply to At <C 1, because to ^> |<7|; (2) the condition |24|) 
contains the square root of n, whereas the eigenvalues 
of exp(±zumAi) depend on the number n itself, which is 
much greater than y/n if n 1. Consequently, even the 
simplest microscopic model gives rise to a QJS, which 
is, strictly speaking, different from the SD jump super- 
operator, coinciding with the former only for the diagonal 
elements \n) (n\ of the density matrix in the Fock basis. 

If condition (f24|) is not satisfied, we propose that the 
QJS can be defined by averaging the transition super- 
operator (|22[) over the interaction time T, because the 
exact instant within (0, T) at which the photodetection 
occurs in each run is unknown, so a reasonable hypoth- 
esis is that these events happen randomly with uniform 
probability distribution: 



jT P=-fJ a dt~(t)p. 



Writing the field density operator as 



m,n— 



p mn \m)(n\, 



we have 



JtP = ^2 PrnnVrrmfn ln \m - 1) (n - 1| , 

m,n— 1 



where 



2A 



T 



fmn = TF / e wt{ - n - m ^ Xt S m {t)S n (t) dt. 



(27) 



(28) 



(29) 



(30) 
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It is natural to suppose that the product AT is big 
enough, so that the photodetection can happen with high 
probability. Mathematically, it means that we assume 
that exp(— AT) C 1, If A <C u> (this is also a natural 
assumption), then the off-diagonal coefficients f mn with 
m 7^ n are very small due to fast oscillations of the in- 
tegrand in Eq. (|2"T)) . so they can be neglected (a rough 
estimation gives for these terms the order of magnitude 
0(\/uj), compared with the diagonal coefficients /„„). 
Consequently, the microscopic model leads to the non- 
linear diagonal QJS of the form 



Jp = 7diag \F(h)apa) F(n)] , 



(31) 



where diag(A) means the diagonal part of the operator 
A in the Fock basis. The function F(n) can be restored 
from the coefficients f nn (apart the constant factor which 
can be included in the coefficient 7) as 



F(n) = y/f, 



n+l,n+l- 



(32) 



Under the condition exp(— AT) <C 1, the upper limit 
of integration in Eq. ([27]) can be extended formally to 
infinity, with exponentially small error. Then, taking 
into account the definition of the function S n (t) (TS")) . 
we arrive at integrals of the form 



dte~ xt x 



sin 



W^ 2 



for x < 1 
for x = 1 



sinh {p£)j p? for x > 1 



which can be calculated exactly (see, e.g., Eqs. 3.893.2 
and 3.541.1 from [HI)- The final result does not depend 
on A or x (and it is the same for either x < 1 or % > 1): 



fnn = (nT)~ 

Thus we obtain the QJS 



(33) 



•hp = It Pnn \n - l)(n - 1| = 7 T diag (E-pE+) , 

n=l 

(34) 

where 7^ = T , and the operators E— and E+ are de- 
fined by Eq. ([7]). Notice that, in principle, 7t is differ- 
ent from "fsD- Moreover, the super-operator (|34j) derived 
from the microscopic model turns out to be different from 
the phenomenological QJS ^ studied in (2(J, HH ■ The 
difference is that Jy has no off-diagonal matrix elements, 
while Je has. We see that the microscopic model con- 
cerned (which can be justified in the case of big number 
of photons in the field mode) predicts that each photo- 
count not only diminishes the number of photons in the 
mode exactly by one, but also destroys off-diagonal el- 
ements, which means the total decoherence of the field 
due to the interaction with MA. 

Note, however, that the formula (f33|) holds under the 
assumption that the upper limit of integration in Eq. 
(|3T))) can be extended to the infinity. But this cannot be 
done if parameter x is very big. Indeed, for x > 1 an d 



AT > 1, the integrand in (30]) at t 

( 



to exp 



-AT 1 



T is proportional 
, so it is not small when 



Vi - 'Vx 2 

n/x 2 <C 1- Calculating the integral in the finite limits 
under the conditions n/x 2 "C 1 and AT 3> 1, we obtain 
the approximate formula 



fnn = (Tn)- 1 {1 - exp [—XTn/ (2 X 2 )] } , 



(35) 



which shows that /„„ does not depend on n if 
XTn/ (2x 2 ) < 1. Thus we see how the QJS (33} can be 
continuously transformed to the SD jump super-operator 
(5|); when the number n changes from big to relatively 
small values. It should be emphasized, nonetheless, that 
the off-diagonal coefficients f mn remain small even in this 
limit. Their magnitude approaches that of the diagonal 
coefficients only in the case of A ~ u>, which does not 
seem to be very physical. 



III. MODEL OF HARMONIC OSCILLATOR 
DETECTOR 

Now let us consider another model, where the role of 
the detector is played by a harmonic oscillator interacting 
with one EM field mode. This is a simplified version of 
the model proposed b y M ollow [29| (for its applications in 
other areas see, e.g., [30] and references therein). In the 
rotating wave approximation (whose validity was studied, 
e.g., in Ref. |3l|) the Hamiltonian is 



H = uj a a)a + uJbb^b + gab^ + g*a'b, 



(36) 



where the mode b assumes the role of the detector and 
the mode a corresponds to the EM field (w& and u a are 
the corresponding frequencies and g is the detector-field 
coupling constant). In the following we shall repeat the 
same procedures we did in the sectionHU The dissipation 
effects due to the macroscopic part of the MA, associated 
to the mode 6, can be taken into account by means of the 
master equation in the form 



dp 



+ i 



HeffP~pHl 



eff 



2Xbpbl 



(37) 



with the effective Hamiltonian 



H eff = H - iXtfb = {uj b - iX) tfb 
J rL0 a a) a + gba) + g*b^a. 



(38) 



The evolution operator U(t) = exp(—iH e fft) for the 
quadratic Hamiltonian (|38[) can be calculated by means 
of several different app roaches (32| . Here we use the alge- 
braic approach 0, [33|, HH, [3f| , since Hamiltonian ([38]) is 
a linear combination of the generators of algebra su(l, 1) 

K+ = 6 f a, iv_ = -ba\ K = (tfb - a f a)/2, 



[K , K ± ]=±K±, [K_, K, 



2K . 
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The evolution operator can be factorized as 

U(t) = e - %mN e A{t)K +e B( - t)Ko e c{t)K - , 

where 

N = (b^b + a) a) /2, Q = uj b + to a - iX. 
The time-dependent coefficients are 

= ig* sin(7?t) c ^ _ igsin(7]t) 



V T (t) 



B(t) = -21nT(i). 



rjT(t) ' 



with 



(39) 



(40) 



(41) 



T (t) = cosfai) + i [oJ ba /(2r))] sin(ryi), (42) 



Since the parameter ZU is big, the off-diagonal coefficients 
fmn with n/m are very small due to the strongly oscil- 
lating factor cxp[iZJt(n — to)]. Consequently they can be 
neglected in the first approximation, and we arrive again 
at the diagonal QJS of the form (|3~l]) . 

We notice that the exact analytical expression for the 
integral in Eq. (|48[) is so complicated (even if to = n), 
that it is difficult to use it. For example, in the limit 
X — > 1 Eq. (|48p can be reduced to the form 



fn 



XT/2 



dyy'(l+yY n -'exp(-2ny). (50) 



Replacing the upper limit by the infinity, we recognize the 
integral representation of the Tricomi confluent hyperge- 
ometric function '3/ (a; c; z) [36| . Thus we have (neglecting 
small corrections of the order of exp(— AT)) 



ui ba = ui b - uj a - iX, rj = (\g\ 2 + uJ 2 a /4) 



1/2 



(43) 



Assuming that the detector is in resonance with the EM 
field's mode one gets u ba = —iX and 

T (i) = cosfot) + [A/(2?7o)] sin(riot), (44) 



Vo = (\g\ 2 A 2 /4) 



1/2 



(45) 



If, initially, the detector oscillator is in the ground state 
|0b), the time-dependent transition operator, correspond- 
ing to the absorption of one photon from the EM field, 
defined in IpOjl . is 



r(t) = (h\u(t)\o b ) 



A(t) exp 



— (iflt + B(t)) (a + o + l) 



a (46) 



and the transition super-operator becomes 

1 



E(t)p = 2A|A(t)|^exp 



(iQt + B(t))(a?a+1) 



x apa) exp 



-(iQ*t- B*(t))(a + a + l) 



.(47) 



For "small" t = At and few photons in the cav- 
ity, the QJS (f2"5|) is recovered. Considering, instead, 
the time-averaged QJS, one has Eqs. i[2T ]) -([29 |) . For 
X = X/(2\g\) < 1 (when the parameter r]o is real) one 
can represent the coefficients f mn as (we consider the 
resonance case with uj a = u> b = u>) 



where 



4* 



dz [cos(z) + £sin(z)] 



m+n—2 



T(i-x 2 y' 2 J Q 

x sin 2 (z) exp [iZJz(n — to) — £z(m + n)] , (48) 



J? I 



T 



*(3;2n + 2;2n). 



(51) 



Although the ^-function in the right-hand side of Eq. 
(|5ip can be rewritten in terms of the associated Laguerre 
polynomials [36[ as 



*(3;2n + 2;2n) 



2(2n) 



(2tQ! r(-l-2n) f g x 



(52) 



neither Eq. (p)T]) nor Eq. (|5"2")l help us to understand the 
behavior of the coefficient f nn as function of n. Therefore 
it is worth trying to find simple approximate formulas for 
the integral in (|48l) . 

If x ^ lj then also ( < 1, so we can neglect the term 
^sin(z) in the integrand of Eq. (|48[) and the function 
sin 2 (2:)[cos(2;)] 2 ™ _2 can be replaced by its average value 
taken over the period 2ir of fast (in the scale determined 
by the characteristic time £ _1 ) oscillations. After simple 
algebra we obtain (replacing the upper limit of integra- 
tion Z by infinity) 



fnn 



4(2n-2)! 



T(2"n!) 2 



X< 1. 



(53) 



Using Stirling's formula nl f=a \/2im(n/e) n , we can write 
for n ^> 1 



fn 



TV im 5 



(54) 



This function corresponds to the QJS (f3lj) with 

F(n) = F 5 (n) = (n+1)- 5 / 4 , 7 = 75 = {T^y 1 . (55) 

Thus, differently from the case of two-level detector, in 
the simplest version of the oscillator detector model the 
lowering operator contains the factor (n + l)~ 5//4 , instead 
of (n + l) -1 / 2 as in the "E- model" ([6]) or simply i as in 
the SD model ©. 

The case x "C 1 is not very realistic from the practical 
point of view, since it corresponds to the detector with 
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very low efficiency. However, we can calculate the inte- 
gral (|48p with arbitrary £ approximately, assuming that 
n 1 and using the method of steepest descent. Rewrit- 
ing the integrand as exp[G(z)], one can easily verify that 
the points of maxima of the function 

G(z) = 2 ln[sin(z)] + 2(n - 1) In [cos (z) + £ sin(z)] - 2£nz 

are given by the formula z k = ±zo + kit, where 

z =tan" 1 (^), fj, = (fn + n - l)~ 1/2 , (56) 

k = 0, 1, 2, . . . for the plus sign and k = 1, 2, . . . for the 
minus sign. One can verify that 



cxp [G(z k )} 



M 2 (l + ^) 2 " 
(1 + /J 2 )™ 



■ exp (— 2zo£n — 2t;imk) 



(57) 

The second derivatives of the function G(z) at the points 
of maxima do not depend on k: 



G"{z k ) 



4n(g 2 + 1) 



(58) 



Using Eqs. (|57p and (|58p and performing summation over 
k we find (taking Z — oo) 

, xVRl + £ M ) 2 "~ 3 / 2 exp (-2z £n) 

(59) 

Although the application of the steepest descent method 
can be justified for n ^> 1, formula (|59"|) seems to be a 
good approximation for n ~ 1, too. For example, for 
n = 1 (when /i = £ _1 ) it yields 

r/n w 4x\/^coth«) cxp [-2£ tan" 1 (£" 1 )] , (60) 

and the numerical values of (|60p in the whole interval 
< % < 1 are not very far from the exact value T/ n = 1, 
which holds independently of x, as far as the upper limit 
of integration in (|48[) can be extended to the infinity. 
For n^> 1 (when jt< 1) Eq. ([59]) can be simplified as 



fnn (x) 



XV8tt 
eT 



~ 3/2 coth 



X < 1. (61) 



For x "C 1 the function (j6"Tj) assumes the form (fS"4"|) , with 
slightly different coefficient 7' = (eT)~ 1 y/8/n m I.O475. 

For x > 1 (when parameter 7/0 is imaginary) we have, 
instead of (|48|) . the integral (considering diagonal coeffi- 
cients only) 



4y f Y 
fnn = j^~2 _ ^3/2 J o [cosh(z) + < sinh(z) 



x sinh 2 (z) exp (— 2n£z) . 



where 



C = X/VX 2 -1, Y = \T/{2(). 



(62) 
(63) 



Applying again the steepest descent method, we have 
now the only point of maximum 



tanh _1 (i/), v = [(C 2 - 1) n+ l] 



-1/2 



(64) 



Taking into account the value of the second derivative of 
the logarithm of integrand at this point, 

4n(C 2 



we obtain 

fnn ' 



i) 



x 



xVMl + <>) 2 "~ 3/2 (l - i/)"^- 1 ) 
T^/Ti{n + x 2 - 1)(1 + ^)™ (C+1) 



1 



(65) 



(66) 



One can check that the limit of formula (|66|) at x — » 1 
coincides with the analogous limit of formula (|59|) . so the 
transition through the point % = 1 is continuous. 

The asymptotical form of (|6"6"|) for n ^> x 2 is the same 
as (|6Tj) . except for the last factor: 



fnn(x) 



XV87 
eT 



-3/2 



X> 1, 



(67) 



Applying the steepest descent method to the integral (|50|) 
(for n 1), we obtain the same result (|6T|) with x = 1. 
Thus for x ~ 1 (not too small and not too big) we obtain 
the QJS in the form ((3TJ> with 



F(h) = F 3 (h) = (n + 1)- 3 / 4 , 7 



73 



eT 



(68) 



For very big values of parameter x (exceeding y/n) the 
steepest descent method cannot be used, because the 
second derivative of the logarithm of integrand, given 
by Eq. (|6"5"|) . becomes small, and because the coordi- 
nate z max , determined by Eq. (|64p , tends to infinity, 
while the upper limit Y of integration in (|62p tends to 
the fixed value AT/2. For x > 1, Eq. (68]) holds for 
the values of n satisfying approximately the inequality 
n > n* ~ 4x 2 exp(— AT). If n < n*, then it can be 
shown that Eq. (|62p leads to the same approximate for- 
mula ([35]) as in the model of two-level detector, so the SD 
super-operator (however, without off-diagonal elements) 
is restored for relatively not very big values of n. 

In Figures [T] and [2] we show the dependence of diagonal 
coefficients /„„ on the number n for different values of 
parameter x> obtained by numerical integration of (|48p 
and (|62p for the fixed value of the parameter AT = 10; in 
figure [3] we compare them with the approximate analyt- 
ical formulas (|59")) and (|6"6"|) . We see that the coincidence 
is rather satisfactory for big values of n, although there 
are some differences for n ~ 1. We also see in Figure [5] 
that the increase of parameter x results in the appear- 
ance of the SD plateau for small values of n, which goes 
into a slope corresponding to the power-law dependence 
for big values of n. The height of plateaus diminishes 
as x~ 2 in accordance with Eq. (f35|) . because big values 
of x correspond (for fixed values of A and T) to small 
coupling coefficient \g\ 2 between the field and MA and, 
consequently, low probability of photocount. 
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ln(n) n 



FIG. 1: Dependence of diagonal coefficients f nn on the num- 
ber n, obtained by numerical integration of (|48|l and (|62p with 
the fixed value AT = 10, for small and moderate values of the 
parameter \ (from below): \ = 0.1, 0.3, 0.5, 0.8, 1.1. 



FIG. 3: Comparison of numerical integration of (|48[) and (|62[) 
with the approximate analytical formulas (|59p and (|66[l for 
X = 0.5, 1.1, 2, 3, 4 (from below). We defined relative error by 

\Jnn Jnn )/ Jnn 




log(n) 



FIG. 2: Dependence of diagonal coefficients /„„ on the num- 
ber n, obtained by numerical integration of (|62[) with the fixed 
value AT = 10, for big values of the parameter \ (from above) : 
X = 5, 10, 20, 40, 70. Notice the appearance of plateaus corre- 
sponding to SD model for initial values of n; for large n they 
are transformed in curves with the slope given by power-law 
dependence. 



IV. CONCLUSIONS 

Here we presented two microscopic models for deduct- 
ing QJS's. In the first one we supposed that the detector 
behaves like a 2-level atom, and in the second - as a 
harmonic oscillator. The main difference between our 
models and previous ones is that we take into account 
the dissipative effects that arise when one couples the ac- 
tual detector to the phototube. This scheme includes the 



spontaneous decay of the detector with originated pho- 
toelectron emission inside the phototube, which is am- 
plified and viewed as macroscopic electric current. Using 
quantum trajectories approach we deduced general time- 
dependent transition super-operator, responsible for tak- 
ing out a single photon from the field. Since it depends 
explicitly on interaction time, we proposed two distinct 
schemes for obtaining time independent QJS's from it. 
In the first case we assumed that the interaction time is 
small and that there are few photons in the cavity; in 
this situation we recovered the QJS proposed by Srini- 
vas and Davies in both detector models. As a second 
scheme, we calculated time-averaged QJS on the time 
interval during which a photon is certainly absorbed; 
as the result, we obtained different non-linear QJS's for 
the 2-level atom model and the model of harmonic os- 
cillator. In particular, we have shown that for quantum 
states with the predominant contribution of Fock com- 
ponents with big values of n, the QJS has the nonlin- 
ear form (|3ip with the power-law asymptotic function 
F{h) = (n+1) - ^. However, the concrete value of the ex- 
ponent (3 is model-dependent. For the 2-level atom model 
we obtained (3 = 1/2, whereas in the model of harmonic 
oscillator the values j3 = 5/4 and (3 = 3/4 were found, 
depending on the ratio between the spontaneous decay 
frequency of the excited state and the effective frequency 
of coupling between the detector and field mode. Also, 
we have demonstrated how the simple Srinivas-Davies 
QJS arises in the case of states with small number of 
photons. Another important result we obtained is that 
the QJS's, when applied to density matrix' non-diagonal 
elements, are null in average in both models due to the 
strong oscillations of the free field terms. 
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